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Abstract 

We investigate the prospects of measuring the strange quark distri- 
bution function at the Tevatron, using W plus charm quark events. 
The W plus charm quark signal produced by strange quark-gluon 
fusion, sg W~c and sg — > W~^c, is approximately 5% of the in- 
clusive -|- 1 jet cross section for jets with a transverse momentum 
PtU) > 10 GeV. We study the sensitivity of the W plus charm quark 
cross section to the parametrization of the strange quark distribution 
function, and evaluate the various background processes. Strategies 
to identify charm quarks in CDF and D0 are briefly described. 



^Talk presented by S. Keller at the workshop "Physics at Current Accelerators and 
the Supercollider", Argonne National Laboratory, June 2 - 5, 1993, to appear in the 
proceedings. 



1 Introduction 



One of the main goals of deep inelastic scattering experiments is to obtain 
reliable measurements of the parton distribution functions of the proton. Re- 
cently, the CTEQIjll] and CCFR|0] collaborations have determined the ratio of 
momentum fractions of strange quark versus u plus d quarks, k = 2S/ (U+D). 
The CCFR collaboration found k ^ 0.4 in its analysis of di-muon events in 
deep inelastic scattering, whereas the CTEQ collaboration with its global fit 
analysis obtained a result of k ~ 1. The difference between the two results is 
at the 3 — 4(7 level. Among current parametrizations of the parton distribu- 
tion functions the MRSDO set0 best represents the CCFR result. The value 
of K, in this set was fixed to k, = 0.5 at = 4 GeV^. As a representative 
of the CTEQ result, we shall use the CTEQIM set in the following. The 
difference between the two parametrizations in k is approximately a factor 
of two at low Q^, and decreases with increasing Q^. This is illustrated in 
Fig. 1, where the ratio of the strange quark distribution functions for the 
CTEQIM and MRSDO set for = 5 GeV^ and = is shown. 



Figure 1: Ratio of the strange quark distribution functions for the CTEQIM 
set and the MRSDO set for two different values of Q^. 
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At large Q^, the strange quarks that originate from gluon sphtting have 
comparable distribution for the two sets because the two gluon distribution 
functions are very similar. 

Here, we point out that the discrepancy between the CTEQ and CCFR 
result could be resolved by a direct, independent measurement of the strange 
quark distribution function. At the Tevatron, such a 

measurement could be carried out by determining the charm content of 

+ 1 jet events. The paper is organized as follows. In Section 2, the signal 
and the leading background processes are studied. In Section 3, different 
experimental techniques to tag a charm quark inside a jet are discussed. Fi- 
nally, in Section 4, we present our conclusions together with some additional 
remarks. 

2 Signal and Background 

Associated W+ charm production proceeds, at lowest order, through sg and 
sg fusion, sg — >■ W~c and sg W^c. The alternative process where the s 
quark in the reaction is replaced by a d quark, is suppressed by the quark 
mixing matrix element Vcd- This suppression is somewhat compensated by 
the larger d quark distribution function, such that the dg — > Wc cross section 
is about 10% of the sg Wc rate. Since the final state is identical for these 
two subprocesses, the sum of the sg and dg contributions will be considered 
as the "signal" in the following. 

The largest background originates from the production of a cc pair in 
the jet recoiling against the W. If only the c, or the c quark, is identified 
in the jet, the event looks like a signal event. Similarly, a bb pair can be 
produced in the jet, and the b or the b quark misidentified as a charm quark. 
Assuming that all the b quarks are misidentified, approximately 75% (20%) 
of the background originates from a cc (bb) pair produced in a jet initiated 
by a gluon. The remaining 5% is due to the production of a cc pair in a 
quark-initiated jet. 

For our subsequent discussion it is convenient to define the following two 
ratios of cross sections: 



and 

signal + background , . 

= W + ljet ' 

where ^^back ground" includes the three background processes mentioned ear- 
her, and 'W + 1 jet" refers to the total inclusive ly + 1 jet cross section. Ri 
represents the ideal situation in which all background events are completely 
eliminated. R2 describes the case where none of the background is removed. 
In practice, the ratio of events with a tag for a charm quark inside of the 
jet over the total number of events will be measured. The result for this 
ratio will fall in between i?i and R2, because specific methods of tagging the 
charm quark inside of the jet will suppress some part of the background, see 
Section 3. There are three possibilities to reduce the background: 

• Charge reconstruction: for the signal, the W and c quark electric 
charges are correlated. For the cc background, the charm quark has 
the wrong charge 50% of the time. Therefore, if the charges of the W 
and of the c quark are determined, the cc background can be reduced by 
a factor of two. Furthermore, events with the wrong charge correlation 
provide a measurement of the background, that could subsequently be 
subtracted. 

• Cut on the charm quark transverse momentum: since more than one 
charm quark is present in the background processes, the average px of 
the charm quark is smaller than in the signal. 

• Flavor identification: if the bottom quark is identified, the bb back- 
ground can be subtracted. 

To numerically simulate the signal and background processes we have used 
the Monte-Carlo program PYTHIA0] (version 5.6). For W + cc and W + bb 
production, we have compared the result of PYTHIA with that of the matrix 
element calculation of Ref. ^. The results of both calculations are in general 
agreement, with PYTHIA resulting in a somewhat larger background. In our 
simulations, a "jet" is defined in the following way. The direction of the sum 
of the momenta of all the partons produced in the shower is taken as the 
center of a cone of radius Ar]'^ + A</)^ = 0.7, where r] is the pseudorapidity 
and the azimuthal angle. All the partons inside the cone are considered to 
be part of the jet. The W is assumed to decay into an electron (or positron) 
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and a neutrino. To simulate the acceptance of a real detector, we impose the 
following cuts on the final state particles: 

Pr(e)>20GeV, |r/(e)| < 1. 

> 20 GeV, (3) 
Pt(j)> 10 GeV, \v{j)\<^- 

Our results for Ri and R2 do not depend sensitively on the cuts. Events with 
two charm quarks inside a jet are counted twice. Numerical values for Ri and 
i?2 are presented in Table |l] for the MRSDO and CTEQIM parametrizations. 



Table 1: The ratios Ri and R2 for the MRSDO and CTEQIM parametriza- 
tion. 





Rii%) 


R2{%) 


MRSDO 


4.2 


8.4 


CTEQIM 


5.3 


9.7 



As can be seen in Table |I|, R2 is about twice as large as Ri , which means 
that the signal-to-background ratio is of order unity. The average values for 
the ratios Ri and R2 are roughly 5% and 9%, respectively. To quantify the 
difference between the two sets, we define the ratio 

_ R{CTEQ1M) - R{MRSDO) 

~ \{R{CTEQIM) + R{MRSDO)) ' ^ ^ 

where R stands for _Ri or R2. A is equal to 23% and 14% for _Ri and R2, re- 
spectively. In order to perform a meaningful measurement, the experimental 
uncertainty in measuring Ri or R2 must be less than the corresponding value 
for A. In Fig. 2, the ratios -Ri and R2 are shown for the two parametrizations 
function of the px of the jet. 
The two sets of parton distribution functions yield the same values for 
the inclusive 1^ + 1 jet cross section and the background to within 1%. The 
variation of the ratios -Ri and R2 with the parton distribution functions 
therefore directly reflects the difference in the strange quark distribution. 
Not surprisingly, both cross section ratios are not very sensitive to changes 
in the factorization scale Q^. Varying between 1/4 and 4 times the default 



5 



Figure 2: The cross section ratios Ri and R2 as a function of the jet pt 
for the MRSDO and CTEQIM parametrizations of the parton distribution 
functions. 

average Q'^ of PYTHIA, the ratios change only by AR/R 4%. The stabihty 
of Ri and R2 with respect to variations in Q'^ indicates that the sensitivity 

of the ratios to the strange quark distribution function is unhkely to be 
overwhelmed by uncertainties originating from higher order QCD corrections. 

Assuming both electron and muon decay channels of the W^, an inte- 
grated luminosity of 10 pb~^ yields about 2000 + 1 jet events for the cuts 
summarized in Eq. (3). This corresponds to approximately 100 W plus charm 
quark signal events, and to about the same number of potential background 
events. From the expected number of signal events it is straightforward to es- 
timate the charm tagging efficiency, e™*", required to be statistically sensitive 
to the variation of the -|- c production cross section with the strange quark 
distribution function. Depending on how efficiently the various background 
processes can be suppressed, wc find e™" ^ 20 — 30% for an integrated 
luminosity of 10 pb~^. Note that e™" scales like (JCdt)^^. 
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3 Charm Quark Tagging in CDF and 



CDF and D0 explore three different strategies to identify charm quarks: 

1. Search for a displaced secondary vertex in the vertex detector. The 
efficiency to tag b quarks with the SVX[^ of CDF is about 10-20%, 
depending on the pt range. The tagging efficiency for charm quarks is 
expected to be smaller than that for bottom quarks, due to the smaller 
mass and decay track multiplicity of the charmed hadrons. 

2. Reconstruction of exclusive nonleptonic charmed baryon or meson de- 
cays. CDF, for example, uses the decay channel Kir to identify 
semileptonic B meson decays^. Other exclusive channels will be added 
in the future, and an efficiency of a few percent should be reached. 

3. Searching for inclusive semileptonic charm decays pj. The average 
inclusive semileptonic charm decay branching ratio is i?(c ^ ev^^v) ~ 
10%. If one assumes a reconstruction efficiency for a muon inside a jet 
of the order of 50%, a total charm tagging efficiency from semileptonic 
charm decays of the order of 5% may well be possible. 

Combined, the three methods may yield an overall charm detection effi- 
ciency of about 10%. 

4 Conclusions 

We have studied the prospects for measuring the strange quark distribution 
function in -|- c production at the Tevatron. The method is similar to the 
one described in Ref. |^ for measuring the charm quark distribution function 
in 7 plus charm production. Our results indicate that, for the data sample 
accumulated in the 1992-93 run, the expected charm tagging efficiencies are 
a limiting factor. For 100 pb^^, however, one could seriously attempt to 
determine the strange quark distribution function from W + c production. 

In our analysis we have concentrated on the charm content of + 1 jet 
events with priJ) > 10 GeV. Alternatively one could search foiW + c pro- 
duction in the inclusive W sample, without requiring the presence of a high 
transverse momentum jet. The advantage here would be a significant in- 
crease of the number of signal events. However, due to the smaller average 
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transverse momentum of the charm quarks in the inclusive W sample, the 
charm quark detection efficiency is expected to be reduced. 

Clearly, more experimental and theoretical work is needed in order to de- 
termine reliably whether a measurement of the s quark distribution function 
in ly + c production is feasible. 

If the strange quark distribution function is measured precisely in other 
experiments, W plus charm quark production may eventually be used to 
measure the quark mixing matrix element Vcs at high and compare it 
with the value extracted from low energy experiments [p^ ] . 
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